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(g) MR-Elastographle-Verfahren 

@ Die Erfindung betrlfft ein MR-Elastographie Verfahren, 
bei dem in einem Untersuchungsobjekt mechanische 
Schwingungen erzeugt werden und die Kernmagnetisie- 
rung angeregt wird. Dabei wird die Auslenkung (U) in ei- 
nem dreidimensionalen Bereich nach Betrag und Phase 
fur drei zueinander senkrechte Richtungen bestimmt und 
wenigstens ein mechanischer Parameter (E) aus diesen 
Werten der Auslenkung und aus deren raumlichen Ablei- 
tungen in wenigstens einem Teil des dreidimensionalen 
Bereichs berechnet. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren zur Bestimmung mechanischer Parameter eines Untersuchungsobjektes mit den 
Schritten: 

5 

a) Erzeugung von mechanischen Schwingungen in dem Untersuchungsobjekt, 

b) Anregung der Kernmagnetisierung in Verbindung mit einem zu den mechanischen Schwingungen synchronen 
magnetischen Gradientenfeld und Empfangen der entstehenden MR-Signale in dem Untersuchungsobjekt zur Er- 
zeugung eines MR-Phasen-Bildes 

10 c) Verandem der Richtung des Gradienten des Gradientenfeldes und/oder der PhasendifFerenz zwischen den me- 

chanischen Schwingungen und dem Gradientenfeld 

d) mehrfache Wiederholung der Schritte a-c) 

e) Bestimmung der durch die mechanischen Schwingungen hervorgerufenen Auslenkung der Kemspins aus ihrer 
Ruhelage anhand der MR-Phasen-Bilder und Berechnung wenigstens eines mechanischen Parameters in Abhangig- 

15 keit von der Auslenkung. 

Solche Verfahren sind als MRE- Verfahren bekannt (MRE = Magnet Resonanz Elastographie). Dabei wird die Tatsa- 
che ausgenutzt, daB die Phase in einem MR-Bild des Untersuchungsobjektes sich infolge der darin wirksamen mechani- 
schen Schwingungen verandert. Das AusmaB dieser Anderungen hangt von der Auslenkung (d. h. der Verschiebung aus 
20 der Ruhelage) infolge der mechanischen Schwingung ab. Somit lal3t sich aus den MR-Phasen-Bildern, d. h. Bildem, die 
die Phase der Kernmagnetisierung darstellen, Information iiber bestimmte mechanische Parameter des Gewebes ablei- 
ten, z. B. iiber die Elastizitat. 

Aus der EP-A 708 340 ist ein solches MR-Elastographie- Verfahren bekannt. Dabei werden zunachst zwei MR-Pha- 
sen-Bilder einer Schicht des Untersuchungsobjekts erzeugt. Der Gradient eines zu den mechanischen Schwingungen 

25 synchronen magnetischen Gradientenfeldes hat bei beiden Bildem die gleiche Richtung, jedoch sind die Phasen dieses 
periodischen Gradienten in bezug auf die mechanische Schwingung um 90° versetzt. Danach werden weitere Paare von 
MR-Phasen-Bildem erzeugt, bei denen der periodische Gradient senkrecht zu der Gradientenrichtung beim ersten Paar 
verlauft. Danach wird die Richtung der mechanischen Schwingung in dem Objekt geandert, und es werden weitere Satze 
von MR-Phasen-Bildem akquiriert. 

30 Aus je einem Paar von MR-Phasen-Bildem laBt sich die Wellenlange fiir die verschiedenen Bildpunkte bestimmen. 
Aus der Wellenlange kann der Elastizitatsmodul (Young's modulus) berechnet werden, wenn die Ausbreitungsgeschwin- 
digkeit der Welle in dem Untersuchungsobjekt und dessen Dichte bekannt sind. 

Bei einem anderen Verfahren, das aus Proceedings of ISMRM 1997, p 1905, Vancouver bekannt ist, wird aus einer 
Folge von MR-Phasen-Bildern die Phase der Auslenkung ermittelt und daraus fiir jeden einzelnen Bildpunkt der Elasti- 

35 zitatsmodul berechnet. 

Den bekannten Verfahren ist gemeinsam, daB sie zufriedenstellende Resultate nur liefem, wenn in dem Untersu- 
chungsobjekt keine Reflexionen auftreten und wenn sich in dem Objekt nur transversale Schwingungen ausbreiten. In ei- 
nem realen Untersuchungsobjekt jedoch, beispielsweise dem Korper eines Patienten, sind Reflexionen unvermeidlich 
und es laBt sich auch keine rein transversale Wellenausbreitung erreichen. Uberdies ist es bekannt, daB longitudinale me- 

40 chanische Schwingungen tiefer in einen Korper eindringen konnen, so daB es an sich erwiinscht ware, daB ein moglichst 
groBer Teil der mechanischen Schwingungsenergie sich longitudinal ausbreitet. 

Aufgabe der vorliegenden Erfindung ist es ein Verfahren der eingangs genannten Art so auszugestalten, daB auch bei 
longitudinaler Wellenausbreitung bzw. bei Reflexionen im Untersuchungsobjekt die mechanischen Parameter des Objek- 
tes zuverlassig ennittelt werden konnen. 

45 Diese Aufgabe wird erfindungsgemaB gelost durch die Bestimmung von Betrag und Phase der Auslenkung in einem 
dreidimensionalen Bereich fiir drei zueinander senkrechte Richtungen und Berechnung wenigstens eines mechanischen 
Parameters aus diesen Werten der Auslenkung und aus deren raumlichen Ableitungen in wenigstens einem Teil des drei- 
dimensionalen Bereiches. 

Die Erfindung bemht auf der Erkenntnis, daB sich die Ausbreitung von mechanischen Wellen in viskos-elastischen 
50 Medien durch eine partielle Differential-Gleichung beschreiben laBt, deren Losung fiir jedes Voxel durch Betrag und 
Phase der Auslenkung in drei zueinander senkrechten Richtungen bestimmt ist, sowie durch die raumlichen Ableitungen 
der Auslenkung. Wenn man diese GroBen fiir jedes Voxel bestimmt hat, kann man die gemessene Werte in die Differen- 
tial-Gleichung einsetzen und daraus mindestens einen der in dieser Gleichung enthaltenen mechanischen Parameter be- 
rechnen. 

55 Es geniigt demnach nicht, fiir lediglich eine einzige Richtung im Raum die Auslenkung zu bestimmen. Ebenso wenig 
reicht es aus, die Auslenkung in alien drei Raumrichtungen nur fiir eine Schicht zu erfassen, selbst wenn man nur fiir 
diese Schicht einen mechanischen Parameter (z. B. den Elastizitatsmodul) ermitteln will. Eine raumliche Ableitung der 
Auslenkung in Richtung senkrecht zur Schicht laBt sich namlich nur bestimmen, wenn die Auslenkung auch in Bereichen 

auBerhalb der Schicht bestimmt wird. Die Auslenkung muB daher in einem dreidimensionalen Bereich bestimmt werden, 
60 d. h. die MR-Phasen-Bilder miissen die raumliche Verteilung der Phase der Kemmagnetisiemng in diesem dreidimensio- 
nelen Bereich wiedergeben. 

Wenn die mechanischen Schwingungen so auf das Untersuchungsobjekt einwirken, daB gemaB Anspmch 2 im we- 
sentlichen longitudinale Schwingungen entstehen, ergibt sich eine groBere Eindringtiefe der Schwingungen, so daB die 
mechanischen Parameter - z. B. des Gewebes im menschlichen Korper - in einem groBeren Bereich bestimmt werden 
65 konnen. 

Die in Anspruch 3 beschriebene Ausgestaltung der Erfindung bewirkt, daB bei jeder Wiederholung der Anregung der 
Kemmagnetisiemng eine genaue zeitliche Zuordnung zwischen den mechanischen Schwingungen und den in Verbin- 
dung mit der Anregung der Kemmagnetisiemng erzeugten Gradientenfeldem erfolgt, insbesondere dem zu den mecha- 
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nischen Schwingungen synchronen Gradientenfeld. 

Die Erfindung setzt voraus, daB eine Reihe von MR-Phasen-Bildern von einem dreidimensionalen Bereich erstellt 
wird, wodurch sich relativ lange MeBzeiten ergeben. Um diese MeBzeiten nicht noch dadurch zu verlangem, daB nach ei- 
ner Anregung das Abklingen der Kemmagnetisierung in dem angeregten Bereich abgewartet wird, werden die Schich- 
ten, aus denen sich der anzuregende dreidimensionale Bereich zusammensetzt, gemaB Anspruch 4 im MultisHce- Verfah- 5 
ren angeregt. 

Eine bevorzugte Ausgestaltung der Erfindung ist in Anspruch 5 angegeben. Es konnen zwar auch andere mechanische 
Parameter berechnet werden, z. B. die Dichte des Gewebes, die Poisson-Zahl oder die Dampfung der Welle durch das 
Gewebe, jedoch ist der Elastizitatsmodul fiir die Diagnose am relevantesten. Die Elastizitat ist der mechanische Parame- 
ter, den ein untersuchender Arzt bei der Palpation des Gewebes erfassen will. 10 

Von besonderem Vorteil ist gemaB Anspruch 6 die Bestimmung des Elastizitatsmoduls bei Mamma-Untersuchungen. 

Anspruch 7 beschreibt eine Anordnung zur Durchfiihrung des erfindungsgemaBen Verfahrens 

GemaB Anspruch 8 haben die Welle und das dazu synchrone magnetische Gradientenfeld einen zeitlich sinusformigen 
Verlauf Es ist zwar auch ein anderer periodischer Verlauf moglich - z. B. sagezahnformig, dreieckformig oder rechtek- 
kig, doch bietet der sinusformige Verlauf demgegeniiber Vorteile. Aus den MeBergebnissen lassen sich bis zu 6 verschie- 15 
dene mechanische Parameter berechnen. 

Die Weiterbildung nach Anspruch 9 gestattet, mehr als 6 verschiedene mechanische Parameter zu berechnen. GemaB 
Anspruch 10 konnen dabei der Elastizitatsmodul und ein weiterer mechanischer Parameter, z. B. der Dampfungskoefh- 
zient berechnet werden, und zwar auch dann, wenn einer dieser beiden Parameter - oder beide Parameter - keine isotrope 
GroBe ist. 20 

Anspruch 11 beschreibt ein fiir das erfindungsgemaBe Verfahren geeignetes Computerprogramm. 

Die Erfindung wird nachstehend anhand der Zeichnungen naher erlautert. Es zeigen 

Fig, 1 ein Blockschaltbild einer Einrichtung zur Durchfiihrung des erfindungsgemaBen \ferfahrens. 

Fig, 2 einen Teil eines dafiir geeigneten Untersuchungsgerats, 

Fig, 3 die zeitliche Lage der mechanischen Schwingungen und der MR-Experimente, 25 

Fig. 4 den zeitlichen Verlauf eines einzelnen MR- Experiments, 

Fig. 5 ein Ablaufdiagramm fiir das MR- Akquisitions- Verfahren, 

Fig. 6 eine schematische Darstellung eines Teils des Auswerteverfahrens und 

Fig. 7 ein Ablaufdiagramm des Auswerteverfahrens. 

In Fig, 1 ist mit 1 ein schematisch dargestellter Hauptfeldmagnet bezeichnet, der in einem nicht naher dargestellten 30 
Untersuchungsbereich ein in z-Richtung verlaufendes stationares und im wesentlichen homogenes Magnetfeld mit einer 
Starke von z. B. 1,5 Tesla erzeugt. Weiterhin ist eine Gradientenspulenanordnung 2 vorgesehen, die drei Spulensysteme 
umfaBt, mit denen in z-Richtung verlauf ende magnetische Gradientenf elder Gx, Gy bzw. Gz mit einem Gradienten in x-, 
y- bzw. z-Richtung erzeugt werden konnen. Die Strome fiir die Gradientenspulenanordnung 2 werden von je einem Gra- 
dienten verstarker 3 geliefert. Ihr zeitlicher Verlauf wird von einem Waveform-Generator 4 vorgegeben, und zwar fiir jede 35 
Richtung gesondert. Der Waveform-Generator 4 wird von einer Rechen- und Steuereinheit 5 gesteuert, die den fiir ein 
bestimmtes Untersuchungsverfahren erf order lie hen zeitlichen Verlauf der magnetischen Gradientenfelder Gx, Gy, Gz be- 
rechnet und in den Waveform-Generator 4 ladt. Bei der MR-Untersuchung werden diese Signale aus dem Waveform-Ge- 
nerator 4 ausgelesen und der Gradientenverstarkeranordnung 3 zugefUhrt, die daraus die fur die Gradientenspulenanord- 
nung 2 erforderlichen Strome erzeugt. 40 

Die Steuereinheit 5 wirkt auBerdem noch mit einer Work-Station 6 zusammen, die mit einem Monitor 7 zur Wieder- 
gabe von MR-Bildem versehen ist. Uber eine Tastatur 8 oder eine interaktive Eingabeeinheit 9 sind Eingaben moglich. 

Die Kemmagnetisierung im Untersuchungsbereich kann durch Hochfrequenzimpulse einer Hochfrequenzspule 10 an- 
geregt werden, die an einen Hochfrequenzverstarker 11 angeschlossen ist, der die Ausgangssignale eines Hochfrequenz- 
senders 12 verstarkt. In dem Hochfrequenzsender 12 werden die (komplexen) Einhiillenden der Hochfrequenzimpulse 45 
mit den von einem Oszillator 13 gelieferten Trager schwingungen moduliert, deren Frequenz der Larmor-Frequenz (bei 
einem Hauptmagnetfeld von 1,5 Tesla ca. 63 MHz) entspricht. Die komplexe Einhiillende wird von der Rechen- und 
Steuereinheit in einen Generator 14 geladen, der mit dem Sender 12 gekoppelt ist. 

Die im Untersuchungsbereich erzeugten MR-Signale werden von einer Empfangsspule 20 aufgenommen und von ei- 
nem Verstarker 21 verstarkt. Das verstarkte MR-Signal wird in einem Quadratur-Demodulator 22 durch zwei um 90° ge- 50 
geneinander versetzte Tragerschwingungen des Oszillators 13 demoduliert, so daB zwei Signale erzeugt werden, die als 
Realteil und als Imaginarteil eines komplexen MR-Signals aufgefaBt werden konnen. Diese Signale werden einem Ana- 
log-Digitalwandler 23 zugefiihrt, der daraus MR-Daten bildet. Aus den MR-Daten werden in einer Bildverarbeitungsein- 
heit 24 MR-Bilder rekonstruiert, die die Kemmagnetisierung im Untersuchungsbereich darstellen. Diese MR-Bilder 
werden nach Anwendung eines Rauschmindemngsfilters auf dem Monitor 7 wiedergegeben und sind Ausgangspunkt fiir 55 
die Berechnung von Bildem (ebenfalls in der Bildverarbeitungseinheit 24), die den raumlichen Verlauf eines mechani- 
schen Parameters, z. B. des Elastizitatsmoduls, darstellen. 

Die Steuereinheit 5 triggert dariiber hinaus einen elektrischen Schwingungserzeuger 25 zur Erzeugung eines sinusfor- 
migen Stromes mit einer festen Frequenz im Bereich zwischen ca. 20 Hz und einigen 100 Hz. Diese Schwingung wird ei- 
nem mechanischen Schwingungserzeuger 26 zugefiihrt und bewirkt, daB ein Kolben mit einer an seinem vorderen Ende 60 
befestigten Platte 27 in eine senkrecht zur Plattenebene periodisch bin- und hergehende, vorzugsweise sinusformige Be- 
wegung versetzt wird. 

Fig. 2 zeigt einen fiir Mamma-Untersuchungen geeigneten Aufsatz auf die Patientenlagerungseinheit 28 des MR-Ge- 
rates, dessen Blockschaltbild in Fig. 1 dargestellt ist. Der Aufsatz umfaBt eine auf nicht naher dargestellte Weise auf die 
Patientenlagerungseinheit abstiitzbare Stiitzplatte 29, auf der sich die Patientin 30 befindet und die mit Ausnehmungen 65 
fiir die Briiste 31 der in Bauchlage befindlichen Patientin versehen ist. Es sind in horizontaler Richtung verschiebbare 
Trager 32 fiir den mechanischen Schwingungserzeuger 26 (nur fiir die linke Seite dargestellt) vorgesehen, mit deren 
Hilfe eine Brust gegen eine mittig angeordnete Begrenzung 33 gedriickt werden kann. Durch Erzeugen einer bin- und 
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hergehenden Bewegung der Kolbenplatte 27 lassen sich dann mechanische Schwingungen in der Brust erzeugen. 

Fig. 3 zeigt die zeitliche Zuordnung zwischen den mechanischen Schwingungen (1. Zeile) und den MR-Experimenten 
(2. Zeile). Mit der Erzeugung der Schwingungen mil der Schwingungsdauer T wird begonnen, bevor durch die MR-Ex- 
perimente MR-Daten akquiriert werden. Dadurch soli erreicht werden, dal3 sich fur die mechanischen Schwingungen ein 
5 stationarer Zustand ergibt, bevor mit der MR-Akquisition begonnen wird. Jedes der aufeinanderfolgenden MR-Experi- 
mente hat eine Dauer NT, wobei N cine ganze Zahl ist. Dadurch wird sichergestellt, daB die mechanische Schwingung im 
bezug auf samtliche MR-Experimente eine definierte Phasenlage hat. 

Jedes der MR-Experimente, von denen eines in Fig. 4 dargestellt ist, umfaBt die Anregung der Kemmagnetisierung 
durch mindestens einen Hochfrequenzimpuls und die Erfassung mindestens eines MR-Signals - in diesem Fall eines 
10 Spin-Echosignals. Die Verwendung von Spin-Echosignalen hat den Vorteil, daB die MR-Signale weitgehend frei von 
Phasenfehlem sind. 

Die 1 . Zeile zeigt die Hochfrequenz- Anregungsimpulse, und zwar zunachst einen 90° Hochfrequenzimpuls RE 1 und 
einen 180° Hochfrequenzimpuls RE 2. Beide sind begleitet von einem Schichtselektions-Gradienten Gz, so daB die Kem- 
magnetisierung jeweils nur in einer Schicht angeregt wird. 

15 Fig, 4 stellt zum Teil den typischen Verlauf eines Spin-Echoexperimentes dar mit einem 90°- bzw. einem 180°-Hoch- 
frequenzimpuls RE 1 bzw. RE 2 (1. Zeile), die von einem Schichtselektions-Gradienten Gz begleitet sind, wobei die 
Phase durch einen Phasenkodierungs-Gradienten Gy kodiert und das entstehende Spin-Echosignal in Verbindung mit ei- 
nem Lese-Gradienten G^ ausgelesen wird. Anstelle einer Spin-Echo- Akquisition konnen aber auch andere MR-Akquis- 
tionsschemen benutzt werden, wie sie z. B. in der EP-A 708 340 beschrieben sind. Die Schichtselektions-, Phasenkodie- 

20 rungs- und Lese-Gradienten mussen nicht mit der z-, y- bzw. x-Richtung zusammenf alien, wie in Fig, 4 dargestellt. We- 
sentlich ist nur, daB sie senkrecht zueinander verlauf en. 

Zusatzlich enthalt das Experiment zwei beiderseits des 180° Hochfrequenzimpulses befindliche, zeitlich sinusformig 
verlaufende Gradienten-Schwingungen Gl und G2, deren Periodendauer der Periodendauer T der mechanischen 
Schwingung entspricht und deren Phasendifferenzen mit der mechanischen Schwingung sich um exakt 180° unterschei- 

25 den. In Fig. 4 verlaufen diese sinusformigen Gradienten Gl und G2 in x-Richtung. Im spateren Verlauf des Verfahrens 
verlaufen sie aber statt dessen in z- oder in y-Richtung. Durch die Gradienten-Schwingungen Gl und G2 wird die Phase 
der durch die mechanische Schwingung ausgelenkten Kernspins geandert, wobei das AusmaB der Anderung von der 
Starke der Auslenkung in Richtung der jeweiligen Gradienten abhangt. 

Fig. 5 erlautert den zeitlichen Ablauf der Akquisition der MR-Daten. Nach der Initialisierung (Block 100) wird der 

30 Kolben 27 in kontinuierliche, sinusformige Schwingungen versetzt, die iiber die gesamte MR-Akquisition anhalten 
(Block 101). Nachdem diese Schwingungen stationar geworden sind, wird zunachst im Block 102 ein erstes MR-Expe- 
riment gemaB Fig. 4 erzeugt und das dabei entstehende MR-Signal (5. Zeile von Fig, 4) erfaBt. Wenn die Dauer NT eines 
MR-Experimentes gering ist, z. B. 50msek, dann ist die Kemmagnetisierung in der zuvor angeregten Schicht bei Beginn 
des nachsten MR- Experiments noch nicht abgeklungen. Deshalb wird das MR-Experiment fiir eine andere - in diesem 

35 Fall in z-Richtung versetzte - Schicht wiederholt (Block 103). Das bedeutet, daB lediglich die Tragerfrequenz der Hoch- 
frequenzimpulse RE 1 und RE 2 so verandert wird, daB die andere Schicht angeregt wird, wahrend die zeitlichen Ver- 
laufe in der zweiten bis fiinften Zeile von Fig, 4 unverandert bleiben. Die Schleife mit den Blocken 102 und 103 wird 
dann so oft durchlaufen, wie Schichten in dem zu untersuchenden dreidimensionalen Bereich vorhanden sind, z. B. 20 
mal, wobei jedesmal eine andere Schicht angeregt wird. 

40 Danach wird der Phasen-Kodiemngsgradient im Schritt 104 geandert, und alle Schichten werden emeut angeregt und 
die darin erzeugten Spin-Echosignale erfaBt. Der Phasenkodiemngs-Gradient wird so oft geandert, wie Phasenkodie- 
mngsschritte erforderlich sind, z. B. 128 oder 256 mal. Danach sind die MR-Daten von samtlichen Schichten akquiriert, 
so daB daraus ein MR-Bild rekonstruiert werden kann, das die Phase der Kemmagnetisierungs-Verteilung in dem durch 
die Schichten gebildeten dreidimensionalen Bereich darstellt. 

45 Nachdem auf diese Weise die MR-Daten fiir das erste MR-Phasen-Bild des dreidimensionalen Bereiches akquiriert 
worden sind, erfolgt im Schritt 105 eine zeitliche Verschiebung zwischen den MR-Experimenten und der mechanischen 
Schwingung und zwar so, daB der Beginn eines MR-Experiments nicht mehr - wie in Fig. 3 dargestellt mit dem NuU- 
Durchgang der mechanischen Schwingung zusammenfallt, sondem demgegeniiber um einen Bmchteil einer Periode T - 
z. B. T/8 - versetzt ist (Schritt 105). Danach wird die Schleife mit den Schritten 102, 103, 104 emeut durchlaufen, so daB 

50 sich ein weiteres MR-Phasen-Bild von dem untersuchten dreidimensionalen Bereich ergibt. Die beiden MR-Phasen-Bil- 
der unterscheiden sich nur in den Bereichen, in denen das Gewebe bzw. die dort angeregten Kernspins in x-Richtung aus- 
gelenkt worden sind. Danach werden weitere MR-Phasen-Bilder erzeugt, wobei der zeitliche Versatz auf 2T/8, 3T/8 . . . 
7T/8 geandert wird. SchlieBlich sind die MR-Daten von acht MR-Phasen-Bildem akquiriert, die in unterschiedlicher 
Weise von der Auslenkung der Kemspins in x-Richtung beeinfluBt sind, die im iibrigen aber identisch sind. 

55 Um auch die Auslenkung in einer anderen Richtung als der x-Richtung erfassen zu konnen, werden die sinusformigen 
Gradienten Gl und G2 statt in x-Richtung gemaB Fig, 4 in einer anderen Richtung erzeugt (Block 106), z. B. in y-Rich- 
tung, so daB - nachdem der beschriebene Durchlauf der Schleifen mit den Schritten 102, 103, 104, 105 vollendet ist - ein 
Satz von MR-Phasen-Bildem des dreidimensionalen Bereiches zur Verfiigung steht, der von der Auslenkung der Kern- 
spins in y-Richtung beeinfluBt ist. SchlieBlich wird dies fiir die z-Richtung wiederholt (die sinusformigen Gradienten Gl 

60 und G2 werden in z-Richtung appliziert), so daB SchlieBlich ein weiterer Satz von acht MR-Phasen-Bildern zur Verfii- 
gung steht, die von der Auslenkung der Kemspins in z-Richtung abhangig sind. 
Danach ist das Akquisitionsverfahren beendet (Block 107). 

Anhand der Fig. 7 wird die Verarbeitung der MR-Daten erlautert. Nach der Initialisiemng im Block 100 werden im 
Schritt 201 - wie schon erlautert - aus den akquirierten MR-Daten MR-Phasen-Bilder erzeugt. Durch die Demodulation 
65 des MR-Signals mit zwei um 90° versetzten Schwingungen Demodulator 22 (Fig, 1) entstehen zwei Signale, die als 
Real- und Imaginarteil eines komplexen Signals aufgefaBt werden konnen. Aus den Realteilen der MR-Signale laBt sich 
ein erstes Bild und aus den Imaginarteilen laBt sich ein zweites Bild rekonstmieren. Die Phase laBt sich fiir jedes \bxel 
aus dem Realteil und dem Imaginarteil errechnen, so daB sich ein MR-Phasen-Bild eigibt. Mit der Berechnung der MR- 
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Phasenbilder kann bereits wahrend der Akquisition der MR-Daten begonnen werden. 

Im Schritt 202 kann die Dampfung der mechanischen Wellen im Untersuchungsbereich ermittelt. werden. Dazu wird 
die Ausbreitung der Welle in dem Gewebe entlang eines Strahls in der Richtung der applizierten mechanischen Schwin- 
gung verfolgt, wobei der Einfachheit halber angenommen ist, daB diese entweder in x-, y- oder z-Richtung auf den Kor- 
per einwirkt. Die Phase eines Punktes auf dem Strahl wird sich dabei infolge der Auslenkung sinusformig andem. Auf ei- 5 
nem anderen, weiter vom Einleitungsort der mechanischen Wellen entfemten Punkt auf dem Strahl, wird sich ebenfalls 
ein sinusformiger Verlauf der MR-Phase einstellen, allerdings wegen der Dampfung der Welle mit verringerter Ampli- 
tude. Aus der Anderung der Amplitude entlang des Strahls kann ein Dampfungskoeflfizient yberechnet werden. 

Im Schritt 203 wird dann fur die einzelnen Voxel des dreidimensionalen Untersuchungsbereichs die Auslenkung nach 
Betrag und Phase fur jede der drei Raumrichtungen ermittelt. 10 

Zur Erlauterung wird auf die schematische Darstellung in Fig. 6 verwiesen. In der oberen Reihe sind drei Satze von 
MR-Bildem schematisch dargestellt, deren Phase von der Auslenkung in x, y und z abhangt. Die Bilder sind vereinfa- 
chend als zweidimensionale Bilder dargestellt, obwohl jedes Bild die MR-Phase in einem dreidimensionalen Bereich 
wiedergibt. In jedem der acht zu einem Satz gehorenden MR-Bilder wird die Phase in dem selben Bildpunkt betrachtet. 
Aus den acht Bildem ergeben sich acht Stiitzstellen fiir den zeitlichen Verlauf der Phase in dem betreffenden Bildpunkt. 15 
Aufgrund der sinusformigen Auslenkung ist der zeitliche Verlauf sinusformig, wie in der zweiten Zeile schematisch an- 
gedeutet. Allerdings kann dem sinusformigen Verlauf jeweils noch ein positiver oder negativer Gleichanteil uberlagert 
sein. 

Es wird dann eine Fourier-Transfonnation liber die durch 8 Stiitzstellen definierte Phase durchgefuhrt, d. h. es wird die 
Sinusfunktion ermittelt, die durch die Stiitzstellen bzw. so dicht wie moglich bei diesen Stiitzstellen verlauft. Diese Si- 20 
nusfunktion hat einen Betrag, der der Auslenkung des Voxels proportional ist und eine Phase, die die zeitliche Lage der 
Sinusschwingung in dem betreffenden Voxel in bezug auf irgendeine Referenzphase definiert. Aus dem Bettag A^ und 
der zeitlichen Phase <p der acht Bilder mit in x-Richtung verlaufenden sinusformigen Gradienten Gl und G2 laBt sich ein 
komplexer Wert Ux fur die Auslenkung in dem betreffenden Voxel wie folgt definieren: 



A^d>=u^=Vx+jWx (1) 



- Y— = — AU + — = V(VU) (2) 
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Vx und Wx stellen dabei den Real- bzw. den Imaginarteil der komplexen GroBe Ux dar. 

Dies wird fiir samtliche Voxel der acht Bilder wiederholt, wonach zwei (fiir einen dreidimensionalen Bereich defi- 
nierte) Bilder entstehen, die den Realteil Vx bzw. den Imaginarteil Wx der Auslenkung darstellen. Das gleiche wird fiir die 30 
Satze von MR-Phasen-Bildem wiederholt, die mit einem in y- bzw. z-Richtung verlaufenden sinusformigen Gradienten 
akquiriert wurden und deren Phase daher durch die Auslenkung in y- bzw. z-Richtung definiert ist. Auf diese Weise er- 
geben sich am Ende des Schrittes 203 drei Bilder, die den Realteil Vx, Vy bzw. Vz der Auslenkung in x-, y- bzw. z-Richtung 
darstellen und drei Bilder, die den Imaginarteil Wx, Wy und Wz der Auslenkung darstellen. 

Der dann folgende Verfalirensschritt 204 beruht auf der Erkenntnis, daB die so gefundenen Werte der Auslenkung Ux, 35 
Uy und Uz Losungen der partiellen Differential-Gleichung darstellen, die die Auslenkung in den betreffenden Voxeln be- 
schreibt. Dementsprechend werden die gefundenen Werte der Auslenkung im Schritt 204 in die Differential-Gleichung 
eingesetzt, um die mechanischen Parameter des Gewebes z. B. seine Elastizitatsmodul berechnen zu konnen. 

Die partielle Differential-Gleichung, die die Ausbreitung einer Welle in einem viskosen, elastischen Medium be- 
schreibt, hat die folgende Form: 40 

d^U du 
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dt^ dt 2(1 + a) 2(l + a)(l-2a) 

In dieser Diffential-Gleichung stellen p die Dichte in dem Untersuchungsbereich, Y die im Schritt 202 ermittelte 
Dampfung der Welle, c die sogenannte Poisson-Zahl (ca. 0,49 fiir menschliches Gewebe) und E den Elastizitatsmodul 
dar. U ist der Vektor der Auslenkung in x-, y- und z-Richtung. V stellt den sogenannten Nabla-Operator dar, d. h. einen 
Vektor fiir den die Beziehung gilt 

A ist der sogenannte Laplace-Operator fiir den bekanntlich die Beziehung gilt: 
A = V • V (4) 
Der An satz 

U = u • ei^t (5) 60 

wobei u ein ausschlieBlich ortsabhangiger Vektor und CO die (Kreis) Frequenz der mechanischen Schwingung darstellt, 
fiihrt nach dem Einsetzen in Gleichung 2 zu der folgenden Diffential-Gleichung: 

u stellt den Vektor der Auslenkung dar, dessen Komponenten in x-, y und z-Richtung sich gemaB der Gleichung 
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U = Ux, Uy, (7) 

zu Ux, Uy und Uz ergeben. Jede Vektor-Komponente laBt sich wiederum als eine komplexe GroBe darstellen, gemaB der 
Beziehung 

5 

Ux=Vx+jWx; Uy=Vy+jWy; Uz=Vz+jWz (8) 

Aus dem Vorstehenden wird deutlich, daB, obwohl Gleichung 6 nur eine einzige DifFerential-Gleichung fiir u darstellt, 
daraus sechs Gleichungen ableitbar sind - je drei fiir jede der drei Richtungen und je zwei fiir Real- bzw. Imaginarteil. 
10 Wenn man annimmt, daB auBer der Kreisfrequenz CO auch die Werte von a, yund p in dem jeweiligen Voxel bekannt sind, 
geniigt es, die gemessenen Auslenkungen in nur eine einzige dieser sechs Gleichungen einzusetzen und daraus den Ela- 
stizitatsmodul E in dem betreffenden Voxel zu berechnen. Die Gleichung 
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-p© v^+YCOw^ = — -( — ^ + — ^ + — ^) +— ( ^+ - + ^) (9) 

2(1 + a) ac^ dy^ dz^' c d^d^ &Ldy d^di 

stellt demgemaB die Gleichung fiir den Realteil der x-Komponente dar. c ist dabei eine Konstante, fiir die die Beziehung 

gilt 

c=2(l+a)(l-2a) (10) 



Aus Gleichung 9 wird deutlich, daB auch in diesem einfachsten Fall die Kenntnis von w^, Vx, Vy, Vz und deren raumH- 
chen Ableitungen in x-, y und z-Richtung erforderlich ist. Es miissen daher auch fiir diesen Fall alle drei Komponenten 
25 der Auslenkung nach Real- und Imaginarteil bzw. Betrag und Phase bekannt sein, ebenso wie die raumliche Ableitung 
dieser GroBen in x-, y- und z-Richtung. Dies setzt voraus, daB die Auslenkungen in einem dreidimensionalen Bereich be- 
kannt sein miissen, selbst wenn man den Elastizitatsmodul nur in einem zweidimensionalen Bereich berechnen woUte. 
Um die raumliche Ableitung bzw. den Differential-Quotienten in x-, y- und z-Richtung ermitteln zu konnen, ist es 
zweckmaBig, den im Schritt 203 ermittelten raumlichen Verlauf der komplexen Auslenkungen v und w (fiir alle drei 
30 Raumrichtungen) durch Hiillkurven zu approximieren, die wenigstens zweimal differenzierbar sind, so daB sich die Ab- 
leitungen numerisch berechnen lassen. 

Durch Einsetzen der auf diese Weise gefundenen Werte fiir die Auslenkungen von Wx, Vx, Vy^ Vz und deren numeri- 
schen Ableitungen sowie den als bekannt vorausgesetzten Werten fiir a, p und y laBt sich mit Hilfe von Gleichung 9 der 
Elastizitatsmodul berechnen und zwar fiir samtliche Voxel des dreidimensionalen Bereichs (bis auf diejenigen, die un- 
35 mittelbar am Rande liegen, so daB die raumlichen Ableitungen fiir diese Voxeln nicht ohne weiteres berechenbar sind). 
Nachdem auf diese Weise im Schritt 204 die Werte E fiir die Voxel des Untersuchungsbereichs berechnet worden sind, 
lassen sich daraus Bilder ableiten, die den raumlichen Verlauf des Elastizitatsmoduls in verschiedenen Schichten des Un- 
tersuchungsbereichs darstellen (Block 205). Danach ist das Verfahren beendet (Block 206). 

Vorstehend wurde davon ausgegangen, daB der Elastizitatsmodul eine isotrope GroBe ist, so daB zur Ermittlung des 
40 Elastizitatsmoduls E lediglich eine der 6 erwahnten Gleichungen zu losen ist. Wenn diese Voraussetzung jedoch nicht er- 
fiillt ist, ist E gemaB 



(11) 



ein Tensor 2. Stufe (d. h. eine Matrix), dessen 9 Tensorelemente Cxx . . . ezz jeweils in das Gleichungssystem einzusetzen 
50 ist. Unter der Annahme, daB dieser Tensor symmetrisch ist, d. h. daB die Bedingungen 

®yx~^y> ®zx~®xzj ^zy~^z (12.) 

erfiillt sind, reduziert sich die Zahl der unbekannten Tensorelemente auf sechs. Diese lassen sich aus den schon erwahn- 
55 ten sechs Gleichungen berechnen, wenn die anderen Parameter wie Poisson-Zahl, Dampfung und Dichte bekannt sind. 

Durch Transformation auf die Hauptachsen entsteht dann ein Tensor, bei dem lediglich die Diagonalelemente von Null 

verschieden sind. AuBerdem ergeben sich aus dieser Transformation die Winkel unter denen die Elastizitat maximal bzw. 

minimal ist. Diese Werte sind ebenfalls von klinischem Interesse. 

Bei dem beschriebenen Rekonstruktionsverfahren wurden gewisse Gewebeparameter als bekannt und konstant im ge- 
60 samten Objekt vorausgesetzt, namlich die Dichte p, die Poisson-Zahl a und der Dampfungkoefhzient y. Die Annahme ei- 

ner konstanten Dichte sowie einer konstanten Poisson-Zahl ist theoretisch gerechtfertigt. Bei dem Dampfungkoeffizien- 

ten hingegen ist es jedoch bekannt, daB er unterschiedlich fiir verschiedene Gewebetypen ist, und daB diese Variationen 

signifikant seien konnen. 

Wenn wenigstens einer dieser Parameter - wie etwa der der Dampfungkoefflzient y- zusatzlich zu den oben erwahn- 
65 ten Tensorelementen des Elastizitatsmoduls ermittelt werden soli, ist dies mit dem bisherigen Verfahren nicht ohne wei- 
teres moglich, well mit den sechs aus den Messungen resultierenden Messungen mindestens sieben Unbekannte ermittelt 
werden miiBten. Um zusatzlich die sechs Tensorelemente des Dampfungkoeffizienten (wobei angenommen ist, daB die- 
ser - ebenso wie der Elastizitatsmodul - lediglich sechs voneinander unabhangige Tensorelemente aufweist) bestimmen 
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zu konnen, wird anstelle einer einzigen sinusformigen Schwingung ein Gemisch von wenigstens zwei sinusformigen 
Schwingungen - sowohl zur mechanischen Anregung des Objektes als auch fiir den Gradienten Gl und G2 - benutzt. 

DemgemaB wird der Oszillator mit zwei Frequenzen gleichzeitig angetrieben, d. h. der zeitliche Verlauf der Oszilla- 
torschwingung S(t) ist 

5 

S(t)=al sin(0)lt)+a2 sin(0)2t) (13) 

wobei al und a2 die Amplituden der Schwingungen mit der Frequenz (Ol und C02 sind. Somit wird auch das Objekt 
gleichzeitig auf zwei Frequenzen angeregt. Da nunmehr im Objekt zwei mechanische Wellen mit unterschiedlicher Fre- 
quenz vorliegen, muB auch der die Bewegung codierende Gradient Gl(t) - vgl. Fig. 4 der die MR-Sequenz bewe- 10 
gungssensitiv macht, aus zwei Frequenzen aufgebaut werden, d. h. 

Gl(t)=bl sin((0lt)+b2 sin(co2t) (14) 

wobei bl und b2 die Amplituden der Schwingungen mit der Frequenz OOl und 0)2 sind. 15 
Mit Vorteil wird dabei 

(02=20)1=20) (15) 

gewahlt. Dadurch kann man die aus den Messungen resultierenden Verlaufe der Phase in den einzelnen Bildpunkten (vgl. 20 
Fig, 6) mit nicht mehr Stiitzstellen (z. B. acht) - bei gleicher Qualitat - bestimmen wie bei einer rein sinusformigen 
Schwingung. Ebenso kann al=a2 und ggf. auch bl und b2 gewahlt werden. 

Das beschriebene Rekonstruktionsverfahren verwendet die mittels MR gemessene OsziUation eines jeden \foxels in- 
nerhalb des betrachteten 3D-Volumens. Diese OsziUation wird tiber den Zeitraum einer Schwingung (T=27C/0) sec) mit 
typischerweise 8 Werten abgetastet. Fine Fouriertransformation dieser Schwingung liefert die gesuchten GroBen, um das 25 
Differentialgleichungssystem zu losen, d. h. Amplitude und Phase bei der Frequenz fur jeden Voxel. Aufgrund der zu- 
satzlichen Anregung bei der Frequenz 20) sowie der Anpassung des zeitlichen Verlaufs des Gradienten sind somit auch 
die vorhandenen Informationen bei der Frequenz 20) verwertbar. Bei der Messung von 8 Werten fiir eine Schwingungs- 
periode 27C/0) liefert die Fouriertransformation Informationen fiir die Frequenzen 0), 2o), 3o) und 4o). 

Weil die mechanische Anregung und die Messung bei zweier Frequenzen erfolgt und die gemischten Informationen 30 
durch die Fouriertransformation nachtragliche getrennt werden, liegen nun die Losungen der Wellengleichung fiir zwei 
Frequenzen vor. Anstelle von Gl (6) erhSlt man somit zwei Gleichungssysteme, 



und 



2(l+a) 2(l + a)(l-2o) ^j^^ 



35 



40 



45 



die miteinander verkettet werden konnen. Dadurch ist es jetzt moglich, nicht nur den Elastiziatstensor, sondem auch den 
Dampfungstensor zu ermitteln. Dabei wird angenommen, daB beide gesuchten GroBen nicht frequenzabhangig sind. 
Dies ist in dem betrachteten Frequenzbereich von ca. 50-400 Hz gerechtfertigt. Es ist aber auch moglich, den Elastiziat- 
stensor mit beiden Gleichungssystemen zweimal zu ermitteln. Durch Mittelung iiber die gefunden Werte wird das Signal/ 50 
Rausch- Verhaltnis in den Bilder verbessert. 

Patentanspriiche 

1. Verfahren zur Bestimmung mechanischer Parameter eines Untersuchungsobjektes mit den Schritten: 55 

a) Erzeugung von mechanischen Schwingungen in dem Untersuchungsobjekt, 

b) Anregung der Kemmagnetisierung in Verbindung mit einem zu den mechanischen Schwingungen synchro- 
nen magnetischen Gradientenfeld (Gl, G2) und Empfangen der entstehenden MR-Signale in dem Untersu- 
chungsobjekt zur Erzeugung eines MR-Phasen-Bildes 

c) Verandem der Richtung des Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 60 
mechanischen Schwingungen und dem Gradientenfeld 

d) mehrfache Wiederholung der Schritte a-c) 

e) Bestimmung der durch die mechanischen Schwingungen hervorgerufenen Auslenkung der Kemspins aus 
ihrer Ruhelage anhand der MR-Phasen-Bilder und Berechnung wenigstens eines mechanischen Parameters in 
Abhangigkeit von der Auslenkung 65 

gekennzeichnet durch die Bestimmung von Betrag und Phase der Auslenkung (U) in einem dreidimensionalen Be- 
reich fiir drei zueinander senkrechte Richtungen und Berechnung wenigstens eines mechanischen Parameters (E) 
aus diesen Werten der Auslenkung und aus deren raumlichen Ableitungen in wenigstens einem Teil des dreidimen- 
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sionalen Bereiches. 

2. Verfahren nach Anspruch 1, gekennzeichnet durch die Anregung von im wesentlichen longitudinalen Schwin- 
gungen in dem Untersuchungsobjekt. 

3. Verfahren nach Anspruch 1, gekennzeichnet durch die mehrfache Wiederholung der Anregung der Kernmagne- 
5 tisierung unter fortlaufender Erzeugung der mechanischen Schwingungen in dem Untersuchungsobjekt, wobei der 

zeitliche Ab stand zwischen zwei Wiederholungen ein ganzzahliges \^elf aches der Periodendauer der mechanischen 
Schwingungen betragt. 

4. Verfahren nach Anspruch 1, gekennzeichnet durch die Anregung von zueinander parallelen Schichten durch ein 
MultisHce- Verfahren, wobei nach der Anregung einer Schicht zunachst andere Schichten angeregt werden, bevor 

10 diese Schicht emeut angeregt wird. 

5. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Berechnung des mechanischen Parameters die Be- 
rechnung des Elastizitatsmoduls umfaBt. 

6. Verfahren nach Anspruch 5, gekennzeichnet durch seine Anwendung fiir Mamma-Untersuchungen. 

7. Anordnung zur Durchfuhrung des Verfahrens nach Anspruch 1 mit einem MR-Gerat, einem mechanischen 
15 Schwingungserzeuger, einer Auswerteeinheit, einem den zeitlichen Verlauf von magnetischen Gradientenfeldem 

bestimmenden Generator und einer das MR-Gerat, den Generator, den Schwingungserzeuger und die Auswerteein- 
heit steuemden Steuereinheit, die so programmiert ist, daB folgende Schritte ausgefiihrt 

a) Erzeugung von mechanischen Schwingungen in dem Untersuchungsobjekt, 

b) Anregung der Kemmagnetisierung in Verbindung mit einem zu den mechanischen Schwingungen synchro- 
20 nen magnetischen Gradientenfeld und Empfangender entstehenden MR-Signale in dem Untersuchungsobjekt 

zur Erzeugung eines MR-Phasen-Bildes 

c) Verandem der Richtung des Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 
mechanischen Schwingungen und dem Gradientenfeld 

d) mehrfache Wiederholung der Schritte a-c) 

25 e) Bestimmung von Betrag und Phase der Auslenkung in einem dreidimensionalen Bereich fiir drei zueinan- 

der senkrechte Richtungen und Berechnung wenigstens eines mechanischen Parameters aus diesen Werten der 
Auslenkung und aus deren raumlichen Ableitungen in wenigstens einem Teil des dreidimensionalen Berei- 
ches. 

8. Anordung nach Anspruch 7, gekenzeichnet durch die Erzeugung einer zeitlich sinusformigen Schwingung durch 
30 den Schwingungserzeuger (26) und den Generator (4). 

9. Anordung nach Anspruch 7, gekenzeichnet durch die Anregung einer ein Gemisch von wenigstens zwei sinus- 
formigen Schwingungen umfassenden Schwingung durch den Schwingungserzeuger (26) und den Generator (4), 
wobei die Frequenzen der beiden sinusformigen Schwingungen vorzugsweise in einem ganzzahligen Verhaltnis zu- 
einander stehen. 

35 10. Anordung nach Anspruch 9, gekenzeichnet durch Mittel (205) zur Berechnung des Elastizitatsmoduls (E) und 

eines weiteren mechanischen Parameters, vorzugsweise des Dampfungskoeffizienten (y). 

11. Computerprogramm fiir eine auf ein MR-Gerat, einen Schwingungserzeuger und eine Auswerteeinheit einwir- 
kende Steuereinheit zur Durchfiihrung des Verfahrens nach Anspruch 1 gemaB folgendem Ablauf: 
a) Erzeugung von mechanischen Schwingungen in dem Untersuchungsobjekt, 
40 b) Anregung der Kemmagnetisierung in Verbindung mit einem zu den mechanischen Schwingungen synchro- 

nen magnetischen Gradientenfeld und Empfangen der entstehenden MR-Signale in dem Untersuchungsobjekt 
zur Erzeugung eines MR-Phasen-Bildes 

c) Verandem der Richtung des Gradienten des Gradientenfeldes und/oder der Phasendifferenz zwischen den 
mechanischen Schwingungen und dem Gradientenfeld 
45 d) mehrfache Wiederholung der Schritte a-c) 

e) Bestimmung von Betrag und Phase der Auslenkung in einem dreidimensionalen Bereich fiir drei zueinan- 
der senkrechte Richtungen und Berechnung wenigstens eines mechanischen Parameters aus diesen Werten der 
Auslenkung und aus deren raumUchen Ableitungen in wenigstens einem Teil des dreidimensionalen Berei- 
ches. 

50 

Hierzu 4 Seite(n) Zeichnungen 
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